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Abstract 

A systematic procedure to derive exact solutions of the associated Lame equation 
for an arbitrary value of the energy is presented. Supersymmetric transformations 
in which the seed solutions have factorization energies inside the gaps are used 
to generate new exactly solvable potentials; some of them exhibit an interesting 
property of periodicity defects. 
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1 Introduction 



The exactly solvable potentials of the one-dimensional Schrodinger equation 
are important due to the fact that the full physical information of the system 
is encoded in a small number of analytical expressions. Moreover, they can 
be used to test the convergence of the numerical methods as well as to be 
the departure point for applying the widely used perturbative techniques. By 
exact solvability we mean that the stationary Schrodinger equation for the 

* Corresponding author. 

Email addresses: davidOf is . cinvestav. mx (David J. Fernandez C), 
gangulyasish@rediffmail.com (Asish Ganguly). 



Preprint submitted to Annals of Physics 



24 May 2006 



involved potential (when it is periodic) admits analytic solutions for energies 
in the allowed bands as well as for energies in the gaps. 

A remarkable fact, the existence of a class of potentials which are intermediate 
between the exactly solvable ones and those just solvable through numerical 
techniques, was realized since the date back to 80s of past century. Nowadays 
they are known as quasi-exactly solvable potentials, and their characteristic 
property is that there exist compact analytical expressions of physical infor- 
mation for only a part of the system spectrum. This means that there are some 
missing features of these models which just can be numerically determined. 
A big effort has been observed for years in identifying the potentials which 
are quasi-exactly solvable, and gradually it started to dominate the conviction 
that this class includes the associated Lame potentials [1-8], among others. 
However, some signs have been recently noticed indicating that the associated 
Lame potentials for some integer values of the parameter pair (m, £) belong as 
well to the exactly solvable class. Recently, we proposed [7] an ansatz through 
which one can implement a fitting procedure to automatically fix the values 
of m, £, the ansatz parameters and consequently the analytic solutions of 
the associated Lame equation. Unfortunately, that procedure was not com- 
pletely systematic in the sense that one has to compute case by case for each 
pair of values of (m, £) to find the general solution by adopting the natural 
modification of the ansatz. 

On the other hand, it seems interesting to enlarge the exactly solvable class of 
potentials departing from a given initial one. Several procedures are available 
to do this, and the simplest one is called supersymmetric quantum mechanics 
(SUSY QM) [9-16]. In this approach there is a differential operator of order k 
intertwining the initial and final Hamiltonians, the main ingredients being k 
seed solutions of the initial Schrodinger equation. These seeds can be chosen 
either physical (as it was done previously when using the ground state) or 
non-physical. When the last one are used, it has been possible to surpass the 
usual restriction of the standard first-order SUSY QM that the new levels will 
be created below the ground state energy of the initial Hamiltonian [16]. 

The supersymmetric transformations started to be implemented just few year 
ago to periodic potentials [1, 17-24]. Similarly as for the non-periodic case, at 
the beginning physical seed solutions were used (band edge eigenfunctions) 
[1,17-19]. However, very soon it was realized that non-physical seed solutions 
could be as well employed. In this way, it was possible to generate either 
periodic potentials (when non-physical Bloch solutions were used) or non- 
periodic ones (with periodicity defects) when general linear combinations of 
the non-physical Bloch solutions were employed [7,20-22]. 

For the associated Lame potentials the supersymmetric transformations were 
applied quite recently [7]. However, the treatment was done for first-order 
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transformations and for particular values of the parameter pair (to, £). It would 
be important to implement the SUSY QM for k > 1 and for general integer 
values of the parameter pair (to, £). 

The main aim of this paper is to present a systematic procedure to derive gen- 
eral solutions of the one-dimensional Schrodinger equation for the associated 
Lame potential with an arbitrary energy E 



which will work, in principle, for any integer values of the parameters to and 
£. Our technique is based on the well known Frobenius method, and we will 
naturally arrive at two separate cases, characterized either by to > £ or by 
m = £. A fundamental conclusion of our treatment is that the associated 
Lame potentials for any integer values of the parameters to and £ are exactly 
solvable. 

In addition, by implementing the supersymmetric transformations (by choos- 
ing the obtained non-physical solutions) we will generate new exactly solvable 
potentials from the associated Lame equation. We will restrict ourselves, by 
simplicity, to first and second-order transformations, but this procedure can 
be continued at will for any order of the intertwining operator. 

Next section will be devoted on the discussion of our procedure at length 
to derive the general solutions of the associated Lame equation (for read- 
ers' convenience a brief introduction about elliptic functions is included in 
the Appendix). We will apply, in the subsequent sections the supersymmetric 
transformations, of first and higher order, to generate new exactly solvable 
potentials, which can have periodic structure or periodicity defects depending 
on how we choose the initial seed Schrodinger solutions. We will end the paper 
with our conclusion. 



2 General solution of associated Lame equation 

Let us make some preliminary remarks about the equation (1). It is defined 
on the full real line (x G R) and the modulus parameter k 2 of Jacobian ellip- 
tic functions belongs to the interval (0, 1). The potential is non-singular and 
periodic of real period 2K or K according to m ^ £ or m = £ respectively, 



where K(k) = ft' 2 d<j>/y]\ - k 2 sin 6. Also it is sufficient to consider m > £ 
for their non-negative integer values (see the discussions in Sec III, Ref. [3]). 



Hij}{x) = [-d 2 x + V(x)\ 4){x) = EiP(x), 
V(x) = m(m + l)fc 2 sn 2 x + £(£ + l)k 2 — - 2 
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The equation (1) may be transformed by applying a coordinate translation 

x-iK' 



x — > z 



, K' = K(k'), k —1 — k , 



(2) 



to Weierstrass form 



dz 2 



+ 



m(m + l)p(z) + 



l(l+l)e 2 e 3 
p{z) - ei 



ij = Eij, i/>(x) = il>(z{x)), (3) 



where 



E = e 3 m(m + 1) + [E - 1(1 + l)]e 3 , e i = e 1 - e u % = 2, 3. (4) 



In equation (3) = a;, uj') is Weierstrass elliptic function of half- 

periods uj = Kj y^, uj' = iK' I ^f&l and the real numbers (ei > e2 > 63) are 
related with Weierstrass invariants through the definition p(uJi) = e^, uji = 
a;, c^2 = oj + a;', cj 3 = uj' . Let us denote two linearly independent solutions of 
equation (3) by ip + (z) , ip~ (z) 1 . Then their product ^(z) = ip + ip~ will be a 
solution of the following third order differential equation 



d 3 ^ 
'dz 1 



m(m + l)p(z) + 



+ l)e 2 e 3 



E 



ei 



-2 



m(m + 1) — 



+ l)e 2 e 3 



p'(z)V = 0, 



(5) 



where, throughout this article the prime will denote differentiation with re- 
spect to its argument. It may be mentioned that in the limit £ = or — 1 
Schrodinger equation (1) reduces to the well-known Lame equation and con- 
sequently the singularity of (5) at z — lo\ disappears, i.e., the singularities 
remain only at the poles of p(z) in the complex z-plane. We will now express 
equation (5) in algebraic form by applying the transformations 



y = 61 p{z \ *(!/) = W) ~ ei]'*- 
e 2 

The equation (5) then reduces to (see the Appendix for more details) 



(6) 



d 3 $ 



d 2 $ 



in which Pi(y) are i-th degree polynomials in y given by 



(7) 



1 For convenience we use this notation to distinguish two solutions and should not 
be confused with the popular notation for intertwined eigenstates. 
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P i {y)=2y 2 {e 2 y 2 -3e iy + e 3 ), 

P 3 (y) = 3y[e 2 (3-2£)y 2 -Q ei (l-£)y +e 3 (l-2£)}, 

P 2 ( y ) = 2{e 2 (3f - m 2 - 6£ - m + 3)y 2 

— [E + ei(9£ 2 -m 2 -9£-m + 3)]y + £(2£ - l)e 3 }, 
P^y) = e 2 (2£ - l)(m + £){m - i + l)y + 2£[E + e x {3£ 2 - m 2 - m)}. 



The differential equation (7) is clearly of Fuschian type having four regular 
singular points at y — 0, l,e 3 /e 2 and oo. It is evident that one can construct 
by Frobenius method, in general, a power series solution around any singular 
point which will be valid in a circle of convergence containing no other singu- 
larity. Such solutions are therefore of local nature. We are interested in global 
solutions, if exist, valid in the whole y-plane. But this means that three local 
solutions around three singular points in the finite part must coincide, which 
is possible iff the Frobenius series terminates after a finite number of terms. 
However, a formal solution around y = can be considered in the form 

oo 

$ = J2a r y r+p . (8) 

r=0 

The indicial equation and the recurrence relations for the coefficients are ob- 
tained in straightforward way as follows 



aofo(p) = 0, ai/ (p + 1) + a /i(p) = 0, (9) 
a r+ 2fo(p + r + 2) + ar+Mp + r + 1) + a r f 2 (p + r) = 0, r = 0, 1, 2 . . . , (10) 

where 



Hp) = e 3 p(p - 1 - 2£)(2p -21- 1), (11) 
h{p) = 2(p - £){ ei [m(m + 1) - 3(p - £) 2 ] - E}, (12) 
HP) = HP -m-£)(p + m-£+l)(2p+l- 21). (13) 

From the recurrence relation (10), it is clear that one can not get a finite series 
corresponding to p = 2£+l and £+1/2. The smallest exponent is, in general, 
non-preferable, since it differs from the greatest exponent by an integer and 
so leads to the solution involving logarithmic terms. But it is known that 
under certain conditions [25] logarithmic terms may not appear in the leading 
solution. We will now investigate this possibility. Note that fo(2£+ 1) = and 
so for r = 2£ — 1, recurrence relation (10) reduces to a two-term relation for 
p = 

a 2e J 1 {2£) + a 2 ^ 1 f 2 {2£-l) = Q (14) 
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and Gt2£+i remains arbitrary up to now. However, the coefficients a r for r < 
2£ + 1 could be determined explicitly in the form 



a r = (-l) r a F r / J] f (s) , r = l,2,...2£, 



s=l 



where F r is an r x r determinant : 



/i(r-l)/ 2 (r-2) 
/o(r-l) fi{r-2) / 2 (r-3) 
/o(r-2) A(r-3) / 2 (r - 4) 







A(0) 



(15) 



(16) 



The first step is to show that F 2 £+\ = 0, which will ensure the absence of 
logarithmic terms from the solution. Close inspection of equations (11-13) 
shows the following interesting properties of fo, fi, f 2 for v — 1, 2, 3, . . . 

h{2i-v) = -fo(v + l), fi(2i-v) = -h(v), / 2 (2€-i/-l) = -f 2 (u). (17) 

It follows at once from equation (16) and (17) that F 2 g + i is skew-symmetric 
and so it vanishes identically. Next, we observe that F 2 £ + i may be expressed 
as F 2f+ i = f 1 (2£)F 2£ — f (2£)f 2 (2£ — l)F 2 £_i, which proves the consistency of 
the constraint relation (14). 

Our final step is to fix a 2 ^+i in such a way that, if possible, the series (8) for 
p = terminates, say, after (N + 1) terms. This is equivalent to impose the 
following additional relations 



a N fi(N) + a N -!f 2 (N - 1) = 0, a N+1 = 0, 



(18) 



f 2 (N) = 0, (19) 

where the integer N is to be determined. The equation (19) clearly gives 
acceptable value of iV for 

N = m + £. (20) 

But then (18) reads 

a m+e f 1 (m + £) + a m+i ^if 2 (m + £ - 1) = 0, a m+e+1 = 0. (21) 

We now consider two cases separately to check consistency of relation (21) 
and, if consistent, to fix a 2e+1 , a 2e+2 , . . . a m+e . 
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case i) m = £ 



In this case relation (21) coincides with (14) and so it is clearly consistent. 
Further equation (21) fixes a 2 e+i = 0. As a result, the series (8) for p = 
terminates after (2£ + 1) terms. 



case ii) m — I + v, v — 1, 2, . . . 

In this case recurrence relations (10) and (21) give v relations for v unknowns 

a 2 i+2fo(2£ + 2) + a 2t+1 h{2£ + 1) + a 2l f 2 {2£) = 0, 
a 2 e+3fo(2£ + 3) + a n+2 h{2£ + 2) + a 2 , +1 f 2 {2£ + 1) = 0, 

a 2t+ „U{2£ + v) + a 2e+ „-ifi(2£ + v - 1) + a 2e+ ^ 2 f 2 (2£ + v - 2) = 0, 
a 2t+ Ji(2£ + + a n+v ^f 2 {2£ + v - 1) = 0. 

Then by back substitution we find 

a 2£+r — a 2e , r — L,Z,...is, [22) 

where D r is the minor of F 2e+u+1 _ r in Laplace expansion of the determinant 
F 2 £ +u+ i. This means that the r x r determinant D r is obtained from F 2 £ +u+ i 
by suppressing (2£ + v+ 1 — r) rows and columns in which F 2 e +u+ i_ r is placed. 
Thus we get a finite series (8) terminating after (2£ + u) terms, since a 2 e +u+ \ = 
a 2 i+v+ 2 = • • • = 0. 

Hence, we have proved that the differential equation (7) possesses a polynomial 
solution of degree (m + £) of the form 

$ = a rV\ a ± 0, (23) 

r=0 

for a special choice of the a r 's. The coefficients a r for r > 2£ may be expressed 
in a compact form, valid for both cases m = £ and m > £ 

a 2e+r = ^(-l) r D m _ e _ r r f[ f 2 (2£ + s)/D m _)j a 2e ,r = 1,2, .. .; D_ r = 0. (24) 

We recall that the a r for r <2£ are given by equation (15), while the rest for 
r > 2£ are to be determined from (24). This means that all the coefficients 
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in the series (23) can be expressed in terms of the normalization constant ao, 
which may be taken as 1. Thus, the product ^(z) of two linearly independent 
solutions ip + (z), ij)~(z) of the associated Lame equation (3) may be written in 
the form (up to some inessential constant factor) 



nr=T W) - p(pr)\ 



(25) 



where p(bi), p(&2), • • • p(b m +£) are the zeros of the polynomial J2T=o a A( e i ~ 
t)/e2] r , a r being determined from (15) and (24). At this moment it is worth 
mentioning that for practical computation of b r one has to invert the tran- 
scendental relation p{b r ) = c r , and so, due to the fact that p(z) is an even 
function, an ambiguity of sign appears in the process, which we shall fix now. 
Let us first express the two linearly independent solutions ip + (z) , ip~ (z) in 
terms of their product ^(z). Note that the Wronskian W(ip + , must be 
non- vanishing and without loss of generality we may set 

W = l=i/> + (il>-)' -(i/> + )'i/>-. 



Dividing the above identity by ^ = ip + ip , 

(ln<T)'-(ln^)' = l 



(26) 



Performing the logarithmic differentiation of the product solution \I> = ^j + ip 
with respect to z, we arrive 



(lnV")'+ (ln^ + )' = (In*)'. 
Thus, adding and subtracting the relations (26) and (27), 



(\nifj~ 



(In*)' 



(ln#)'- 



From (28) it follows readily 

, / \ r z c 

ip ± (z) = ^(z)exp\^j y 



dr 



(27) 



(28) 



(29) 



To fix the sign of b r , we will now differentiate the second relation in (28): 



Substitution for (^ + )'/^ + = (ln^+)' from (28) 



+ 



(30) 
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Now multiplying equation (30) by 4\I/ 2 , and using the associated Lame equa- 
tion (3) for the solution we obtain 



2^>^>" + 1 — ^>' 2 = 4\I> 2 



m(m + l)pu) + . . £ 

p{z) - ei 



(31) 



Noting that 6 r are the zeros of ^(z) [see equation (25)], for the values z = 
b r , r — 1, 2, . . . m + £, equation (31) gives 



z=b r 



= 1. 



(32) 



The ambiguity of signs in the values of b r may now be fixed by selecting the 
convention 

*'\ z = br = +1, (33) 

which can also be expressed in terms of p-functions by evaluating at z — bj 
from (25) as 



z=bj 



[ s){b : ) -r 1 Y II [P(y - = +1, J = l,2,...m + £. (34) 



The remaining job now is to express the integrand in (29) in a form suitable 
for performing the integration. Let us express 1/\I/, with \1> given by (25), as 
a sum of partial fractions in the form 



1 



m+i 



- p(6 r )" 



Then, equating each term from both sides of the above identity for z = 
bi, 62, • • ■ b m+ e, and using the relation (34), it is straightforward to obtain 



1 p'{b r ) 



(35) 



To proceed further we introduce two quasi-periodic functions ((z) and a(z) 
by the definitions 

C(z) = -p(z), [\na(z)]' = ((z), (36) 

with the properties 

C(z + 2ui) = ((z) + 2C(^), ((-z) = -C(z), (37) 



a(z + 2ui) = - exp[2((uji)(z + Ui)]a(z), cr(-z) = -cr(z). (38) 
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These are known as Weierstrass zeta and sigma functions respectively. The 
addition formulae for them are 



> / i \ >/ x . >/ x . i p'(z) - p'{y) 

C{z + y) = C{z) + C{y) + 2p( Z )-p(y)> 



(39) 



a(z + y)a(z - y) = -a 2 (z)a 2 (y)[p(z) - p(y)\. (40) 
For our purpose it is useful to rewrite the addition formula (39) in the form 



C(z-y)-C(z + y) 



- p{y) 



-2C(y), 



(41) 



where we have used the fact that p(z) is an even function, but p'{z), ((z) 
and a(z) are odd. It is now not very difficult to rewrite the equation (35) by 
replacing y in (41) by b r : 



i m+e 

= E IC(* " ~ C(Z + K) + 2C(6r)]. 



(42) 



The advantage of writing 1/^ in the form (42) is that one can express the 
argument of the exponential term in the solution (29) in a compact form by 
exploiting the definitions (36) 



f z dr m+e 



ln^M + 2, C (M 
a\z + b r ) 



(43) 



The factor ^J^f(z) in the solution (29) may also be expressed in terms of sigma 
and zeta functions by using (25) and (40) as (up to some factor) 



exp[£zCM] ^ 



a m (z)a e (z + u 1 ] 



II y a ( z + b r )cr(z - b r 



(44) 



r=l 



It may be mentioned that in the above equation we have eliminated the term 
<j(z — uj\) by using the quasi-periodic property (38). Our solutions of the asso- 
ciated Lame equation (1) for any integral values of m, £ will now be presented 
in their final form 



m+l _/ x—iK' 



±b r ) 



X 



m+l 

€CW =F E C{b r ) 

r=l 



(45) 



It is instructive to consider the limit t — > 0, which will reduce the associated 
Lame equation (1) to the ordinary Lame equation 



[— d x + m{m + l)k sn x]ip(x) = Eip(x). 



(46) 
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The solution of (46) can be obtained from (45) in the limit t — > 



IX 



a 



x~iK' 



r m( x-iK' \ 



exp 



V^3 



(47) 



which is exactly identical with the old result obtained for Lame equation in 
Ref. [26]. Further, it may also be noted that the general solution (45) perfectly 
coincides with our previous results [7] for (m, £) = (1, 1) and (2, 1). In the next 
section we will apply the supersymmetric transformations using generalized 
superpotentials based on the general solution (45). 



3 Supersymmetric transformations 



A simple technique to generate new Hamiltonians H with known spectra from 
a given initial one H is the so-called supersymmetric quantum mechanics 
(SUSY QM) [9-16]. In this procedure the non-null action of a finite-order 
differential intertwining operator B such that 

HB = BE (48) 

onto the eigenf unctions of H provides those of H (which in particular is valid 
for the physical eigenfunctions). In the modern approach to the subject the 
intertwining operator B is of fc-th order, and its simplest variant involves k 
seed Schrodinger solutions of H associated to k different factorization energies, 
namely, 

Hui = 6iUi, (49) 
which are annihilated as well by B, i.e., 

B^ = 0, i = l,...,k. (50) 

The new potential V(x) differs from the initial one V(x) by a term involving 
the Wronskian W(u 1: . . . , u k ) of the k seed solutions in the way 

V{x) = V{x) - 2[ln W(u u u k )f. (51) 

It is worth to notice that this generalized approach to SUSY QM allows to 
surpass the typical restriction of the standard first-order version that the new 
levels have to be created below the ground state energy E of H . This is 
possible since the SUSY transformations can be implemented using either 
physical Schrodinger solutions (as it was typically done, e.g., through the use 
of the ground state) or non-physical ones [7] . 

The SUSY QM was applied for a long time to Hamiltonians with a discrete 
part of the spectrum (see e.g. the collection of articles in [10]). Except by 
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few previous works [27], however, the SUSY transformations started to be 
implemented in a systematic way just recently to potentials which are peri- 
odic, specifically to the Lame potentials of (46) [1, 17-24]. Similarly as for 
the non-periodic case, for the Lame equation (46) the SUSY transformations 
were implemented first by using seed physical solutions (the band edge eigen- 
functions of H) [1,17-19]. However, very soon it was understood that more 
general eigenfunctions of H (non-physical ones included) provide once again a 
generalized scheme [20-22]. In particular, it includes the possibility of gener- 
ating either periodic or non-periodic (with periodicity defects) SUSY partner 
potentials, which depends on the use of Bloch solutions (47) or their general 
linear combinations for a fixed set of factorization energies. 

Concerning the associated Lame potentials (1), the SUSY QM has been ap- 
plied recently [7] . The treatment was restricted to first-order transformations, 
the parameter pair (m, £) taking the values (1, 1) and (2, 1). In this paper we 
will implement the supersymmetric transformations to the associated Lame 
potentials for any integer value of the pair (m,£). For the sake of simplicity, 
we will restrict ourselves to first and second-order transformations. 



3.1 First- order SUSY QM 



Let us apply then the SUSY QM of first order to the associated Lame po- 
tentials. We denote by u(x) the involved seed Schrodinger solution and by e 
the corresponding factorization energy. In order to avoid the singularities in 
V(x) we will suppose that e < E . This means that we need to know the exact 
information of the ground state in order to choose e appropriately. Fortunately 
this information is already available in the literature [1-3]. 

We will use in the first place any of the two Bloch solutions (45) as transfor- 
mations functions, i,e., we take u(x) = ^(x). Therefore: 



r, / \ i // m . £(2e 2 + e 2 e 3 ) sn 2 (x — iK') 

[\nu(x)}" = — 7 77J r+£ + 



sn 2 (x — iK') e\ cn 2 (x — iK') 



m+l -y 



§ Sn 2 {X ± yfeb r - iK')' ^ 



By employing that k 2 sn 2 (x — iK') = l/sn 2 a; and k 2 cn 2 (x — iK') = — dn x/sn 2 x 
one arrives to: 

[In u(x)\" = mk 2 sn 2 x + £ - £ ( 2e i + ^ _ y sn V x ± (53) 

e|dn x ^ 

Thus, the periodic first-order SUSY partner potentials generated through the 
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Bloch solutions (45) become: 

2 m-\-t 

V±(x) = m(m - l)k 2 sn 2 x + £(£ - l)k 2 ^- + 2k 2 V sn 2 (a; ± v^A)- (54) 

an x fr{ 

Since the intertwining operator B transforms bounded Schrodinger solutions 
into bounded ones, unbounded into unbounded, etc., it turns out that the 
new potentials have exactly the same band spectrum as the associated Lame 
potentials. 

On the other hand, if we use as seed a nodeless linear combination of the two 
Bloch solutions (45) associated to a given factorization energy e < E , 

u(x) = Aip + (x) + Bip-(x) = Aip + (x)(f) + (x) = Bip-(x)(j)-(x), (55) 

where A + = B/A, A_ = A/B, and 

<f> ± ( x ) = 1 + A ± — — = 1 + A± 11 K , exp - ir - , (56) 

it turns out that the SUSY partners for the associated Lame potentials become 
now: 

V np (x) = V ± (x)-2[ln<f> ± (x)\", (57) 
where V±(x) are given by (54). 

Let us notice that V np (x) is non-periodic in the full real line, a behavior 
characterized specifically by the term [In (f>±{x)]" . In fact, there is a periodicity 
defect in a finite region of x, but the potential acquires an asymptotic periodic 
behavior when we move far away of that domain. Let us remark that in the 
direct approach, even if we would know how to do it, it would be difficult 
to determine the spectrum of the Hamiltonians associated to the non-periodic 
potentials V np (x). However, due to its periodic behavior for large \x\ and since 
the intertwining operator maps bounded Schrodinger solutions into bounded 
ones, unbounded into unbounded, etc., it turns out that the SUSY approach 
provides in a simple way the spectra of the new Hamiltonians, which contain 
the allowed energy bands of the initial associated Lame potential. Moreover, it 
can be shown that l/u(x) is square-integrable, meaning that V np (x) acquires 
an extra isolated bound state at e. 

3.2 Second-order SUSY QM 

Let us apply now the second-order SUSY QM by using two Bloch solutions 
of the form (45) associated to two different factorization energies e±, e 2 , where 
we denote by b r and b' r the corresponding constants. In order to avoid the 
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zeros in the Wronskian, which would produce singularities in V(x), we choose 
e±, €2 to be in the same forbidden gap of H (for a concrete proof see [20]); this 
includes of course the possibility for both to be in the infinity gap below E 
but it allows as well to use solutions in finite gaps (above E ). Once again, this 
means that we will need information about the positions of the band edges, 
which will be taken from earlier works [1-3]. We will select, for definiteness 
both solutions with the upper signs; any other signs combination will lead 
essentially to the same expressions which will be derived below. Thus we take 



u 2 (x) = ip£(x) 



'x-iK' 



cxp 



m+l 



{4=[*CM-£C(6r)]} 

* 63 



r=l 
m+£ 



^{-^[%)-5;^)]}(58) 



r=l 



An appropriate expression for the Wronskian reads: 



W(u 1 ,u 2 ) = g(x) u x (x) u 2 (x), g(x) 



In 



(59) 



where, by using the previous formulae for Ui(x), U2(x), it can be shown that: 



9{X) = -j= E M^TMT + b 'r) ~ W) " + b r) + C(6r)] (60) 

V 63 r= \ 



'e 3 ' x V e 3 

Thus, the modification to the potential is given by: 



[In W( Ul , u 2 )}" = (In g)" + (m^)" + {\nu 2 )" 

= (Jug)" + 2mkhn 2 x + 21- 2 £ ( 2e j+ ^ 



e%dn 2 x 



m+e 



-k 2 E [sn 2 (x + \fe^b r ) + sn 2 (x + \[e$r 



r=l 



(61) 



where we have used (53) to simplify this expression. Finally, the periodic 
second-order SUSY partner potential V(x) of V(x) is given by: 



2 

V(x) = m(m - 3)kVx + £(£ - 3)k 2 ^- 

dn x 

m+l 



+2k 2 J2 [sn 2 (x + yfeb r ) + sn 2 (:r + y/E^V r )] - 2 (In g)" . (62) 



r=l 



On the other hand, if we use as seeds two general linear combinations of the 
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Bloch solutions associated to ei, e 2 , which up to unimportant constant factors 
can be expressed as 



M x ) = *l>i + = ^1 <H, M x ) = i>t + A 2^2 = ^2 <f>2 > 



(63) 



where 



>+ = 1 + Ai 



i>2 



m+e <r(^|p - b r ) 



exp 




(64) 



we will arrive at non-periodic second-order SUSY partner potentials. Indeed, 
a convenient expression for the Wronskian reads: 



W( Ul , u 2 ) = Ul u 2 g np = ^2^2 9 np , 



(65) 



where 



np 




+ 




In ^ 



9 + 




In ^ 



(66) 



with g given by (60). The modification of the potential is thus given by: 



[\nW( Ul ,u 2 )}" = [ln(#+^ 2 + ) 



+ 




(67) 



With this equation it is simple to see that the non-periodic second-order SUSY 
partner of the associated Lame potential becomes: 



V np (x) = V(x) - 2 




(68) 



where V(x) is given by (62). Let us notice that this part of V np (x) coincides 
with the periodic second-order SUSY partner potential previously derived. 
The non-periodicity in V np (x) is produced by the second term in the RHS 
of equation (68), but now this term is more involved than in the first-order 
case. However, the general behavior of V np (x) is quite similar, i.e., there is 
a finite region in the x-domain in which a periodicity defect is produced. 
Hence, similar to the case for first-order SUSY we have obtained a potential, 
which is asymptotically periodic. It is worth to mention that the spectrum of 
V np (x) contains again the allowed energy bands of the initial associated Lame 
potential, but now there will be two extra bound states at e±, e 2 which could 
be interesting for physical applications. 
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Fig. 1. Periodic first-order SUSY partner potential V+(x) (black curve) isospectral 
to the associated Lame potential (gray curve) with m = 3, i = 1, k 2 = 0.95, 
generated by using the Bloch solution u(x) = ip + (x) with e = 4.75. 



4 Applications 



In a previous paper [7] we have found solutions of the stationary Schrddinger 
equation with an arbitrary energy (either physical or non-physical) for the 
associated Lame potentials, the pair (m,£) taking the values (1, 1) and (2, 1). 
We applied there also the first-order SUSY techniques in order to generate new 
potentials with known spectra. Here, we will illustrate our previous general 
results with a different case characterized by (m, £) = (3, 1). For this associated 
Lame potential explicit expressions for some band-edge eigenf unctions and 
eigenvalues are known [3] 2 : 



■0(1) = cnx dn 2 x, 
ip(2) = sna; dn 2 a;, 



E 1 
E2 



1 + 4k 2 , 
1 + 9k 2 , 



^(3,8) 



snicni 
dnx 



sn 2 x - — ( k 2 + 3 ± Vk 4 - 9k' 2 
ok 2 v 



£3,8 = 10 + 2k 2 =f 2Vk 4 - 9k 12 . 
Three other band edge eigenfunctions can be given: 



_ sn 4 x {9k 2 + 16 - Ei) sn 2 x Ef — 2(5k 2 + 18)^ + 9k A + 156/c 2 + 320 
^ dnx 10k 2 dnx 15/c 4 dnx 

2 We will denote here by the band edge eigenfunctions to distinguish them from 
the solutions used in the SUSY transformations. Some mistakes in the ordering of 
levels in Ref. [3] have been corrected here. 



16 



where the eigenvalues E^ i — 0,4,7 are the ordered roots of the cubic equa- 
tion: 

E 3 - (Ilk 2 + 20)E 2 + (19k 4 + 21Qk 2 + QA)E - (9k 6 + 388A; 4 + 448A; 2 ) = 



10 




I I 



-8K 



-4K 







4K 



8K 



Fig. 2. Non-periodic first-order SUSY partner V np (x) (black curve) of the associated 
Lame potential (gray curve) for m = 3, t = 1, k 2 = 0.95, which was generated by 
using u(x) = ip + (x) + ip~ (x) with e = 4.75. The new Hamiltonian H np has an extra 
bound state precisely at e = 4.75. 

Our first task is to evaluate four constants a±, a,2, 03, 04 (note that m + £ = 4) 
in (23), where without loss of generality, we choose oq — 1. Let us first write 
down the basic elements /j for m — 3, £ — 1 from (11-13) 



f (p) = e 3 p(p - 3)(2p - 3), Mp) = 2(p - l){3ei[4 - (p - l) 2 ] - E}, (69) 
/ 2 (p) = e 2 (p - 4)(p + 3)(2p - 1), E = e 3 (E - 2) + 12e 3 . (70) 



Note that the two constants a 1; a 2 have to be determined from (15) 



/ (1)' a2 / (l)/o(2)' 
while the remaining two are to be computed from (22) 



03 = 



£1/2(2) 



-a 2 , 



0,4 = 



£o / 2 (2)/ 2 (3) 
Do 



(71) 



(72) 



Now from (16), one may write quite straightforwardly the Fj's 



^i = /i(0), 



/i(l) / 2 (o) 

/o(l) /i(0) 



(73) 
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and to obtain the ZVs we need 



/i(4) / 2 (3) 
/o(4) /i(3) / 2 (2) 
/ (3) A(2) / 2 (1) 
/ (2) ^(l) / 2 (0) 
/ (1) /i(0) 

Then the final quantities are 

D Q = [minor of F 5 in F 5 ] — 1, 



(74) 



(75) 



£>i = [minor of F 4 in F 5 ] = /i(4), 



-D 2 = [minor of _F 3 in F 5 ] 



/i(4) / 2 (3) 
/o(4) /i(3) 



Finally, using the ingredients (69-77) we obtain 



(76) 
(77) 



a 



9ei — E 



e 3 



a 2 = 



Q3 = 



45(£ + 15ei)el 



(E 2 + 15 ei E + 25e 2 e 3 )e 3 ' 



04 



6e 2 

e 3 

225 4 



{E 2 + 15 ei E + 25e 2 e 3 )e 3 



(78) 



We employ these coefficients then to find the roots c r ,r = 1,...,4 of the 
fourth-order equation 

' e\ — V 



4 

r=0 



e 2 



= 0, 



(79) 



which can be analytically determined, but their explicit expression is too in- 
volved to be shown here. These roots are used then to invert the transcendental 
equation p(b r ) = c r to determine the 6 r 's (with the restriction ^'1^=6,, > 0), 
which are thus inserted in the explicit expressions for ^(x). Finally, the re- 
sulting Bloch solutions can be used, either directly or in the corresponding 
Wronskian, to derive the periodic SUSY partner potentials V±(x) of (54) or 
V(x) of (62). On the other hand, different linear combinations of kind (55) 
or (63) can be used to derive the potentials V np (x) of (57) or (68) which 
have periodicity defects. The final results of these procedures are illustrated 
in figures below. In the four figures we show in gray the original associated 
Lame potential for m = 3, £ = 1, k 2 = 0.95. In figure 1 we show as well in 
black one of its periodic first-order SUSY partners generated through a Bloch 
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solution with e = 4.75 while in Figure 2 it is illustrated one of its non-periodic 
partners for the same e, which is lesser than but close to the lowest band edge 
Eq = 4.79991. On the other hand, in figures 3 and 4 we have drawn similar 
graphs (black curves) for the corresponding second-order SUSY partners, pe- 
riodic and non-periodic respectively. For the periodic case (Figure 3) we have 
used two Bloch solutions ui(x) = ipf(%), u 2 {x) = i>£( x ) associated to the pair 
of factorization energies e\ = 9.4, and €2 = 9.5 which are in the first finite 
energy gap (4.8,9.55). For the non-periodic case we have used the same pair 
of factorization energies, with linear combinations Ui(x) = ipf(x) +ipi(x) and 
u 2 (x) = 1P2 ( x ) ~ 2^2" ( x )- 111 both non-periodic cases the periodicity defects 
are clearly detected. 



10 



5 







j 




1 


1 


1 




If] 



















-4K 4K 



Fig. 3. Periodic second-order SUSY partner potential V{x) (black curve) isospectral 
to the associated Lame potential (gray curve) with m = 3, i = 1, k 2 = 0.95, 
generated by using e± = 9.4, u±(x) = ip^(x) and e 2 = 9.5, U2{x) = ^(x). 

5 Conclusions 

In this article we have shown, by finding explicitly the two solutions of Bloch 
type associated to the stationary Schrodinger equation for an arbitrary value 
of the energy, that the associated Lame potentials for any integer values of the 
parameter pairs (m, £) are exactly solvable. This point is clarifying because 
in most of the works concerning Schrodinger equation with periodic poten- 
tials typically one looks for just the band edge eigenf unctions, thus inducing 
the idea that these solutions are the only ones which could be analytically 
determined. The solved problem is very important for implementing the su- 
persymmetric transformations, of first or higher order, since the non-physical 
Schrodinger solutions can be used as seeds to generate new exactly solvable 
potentials. Consequently, we have generated in this way potentials which are 
either strictly isospectral to the initial one or with some isolated bound states 
embedded at the energy gaps of the initial Hamiltonian. The arising of these 
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Fig. 4. Non-periodic second-order SUSY partner V np (x) (black curve) of the associ- 
ated Lame potential (gray curve) with m = 3, £ = 1, k 2 = 0.95, which was generated 
by using e\ = 9.4, u\{x) = ipf(x) + ipi{x) and 62 = 9.5, U2{x) = 1^2 ( x ) ~ ^2 ( x )- 
The Hamiltonian H np has the same band spectrum as the initial associated Lame 
potential plus two isolated bound states at ei and €2- 

two different kinds of spectra for the new potentials depends on either we 
choose as seeds directly the Bloch solutions or general linear combinations for 
the chosen factorization energies. The potentials with bound states embedded 
into the finite gaps may be interesting for physical applications since the cor- 
responding levels could work as auxiliary transition energies for the electron 
to jump easier between allowed bands. 
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A Appendix 



In the following we will provide a short introduction about elliptic functions 
(for more details, see [28,29]). Three Jacobian elliptic functions are defined by 

sn(x, k) = sin <p, cn(x, k) — cos <p, dn(x, k) = dip/dx, (A.l) 
where amplitude function <p(z, k) is defined by the integral 



i((p,k) = / 
Jo 



dr 



\J\ — k 2 sin 2 r 



(A.2) 



The square of the real number k is called elliptic modulus parameter and 
k 2 e (0,1). k' 2 = 1 — k 2 is called complementary modulus parameter. For 
simplicity, in the text we suppress the explicit modular dependence and write 
snx, cnx, dnx, etc. These are doubly periodic functions of periods 4K, 2iK'; 
4K, AiK' and 2K, AiK' respectively, where the quarter-periods K and K' are 
the real numbers given by 



K(k) = K = z(7r/2,k), K'(k) = K' — K(k'). 



(A.3) 



K is called complete elliptic integral of second kind. Some relevant relations 
are 

sn(x + K) = cn(x + K) = -k'^, dn(x + K) = (A.4) 



dnx ' 



dnx ' 



dnx' 



sn(x + 2K) = -snx, cn(x + 2K) = -cnx, dn(x + 2K) = dnx, (A. 5) 



x+iK') = , cn(x+tK') = , dn(x+tK') = -% , (A. 6) 

/csn x k snx snx 



sn x + cn x = 1, dn x + k sn x — 1, 
and the rules of differentiation are 



(A.7) 



sn'x = cnx dnx, cn'x = — snx dnx, dn'x = — k 2 snx cnx. (A. 8) 



Weierstrass elliptic function p(z; g<i, g-z) = p(z) is defined by 



m,n 



(z — 2muj — 2nu') 2 (2mw + 2nu'Y 



(A.9) 
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where the symbol J2' means summation over all integer values of m, n except 
m = n = 0; u, u' being half-periods of p(z). The invariants g 2 , # 3 are given 
by 

92 = 60 ]T' ? 93 = 140 E' 7 7^6- ( A - 10 ) 

The three numbers ej, % = 1,2,3 are defined by p(u>i) = e^, where iO\ = ou, u; 3 = 
u/ and w 2 = w+o/. To get the relation between Jacobian elliptic functions and 
Weierstrass elliptic function, it is necessary to define uj,uj' in terms of K,K'. 
We have taken the following definition 

cu = — =, u' = -=, A.ll 

which corresponds to the case when the discriminant A = p| — 27g 3 > 0. 
This means that the numbers e$ are always real and consequently they can be 
ordered as ei > e 2 > e 3 , because these are roots of the equation 

At 3 - g 2 t - g 3 = 0. (A. 12) 



The relations between p{z) and snz, cnz, dnz may then be written as 

( z \ _ cn 2 z _ dn 2 z 1 

P -7= = ei + e 3 ^- = e 2 + e 3 ^- = e 3 + e 3 — =-. (A.13) 
V \/e 3 / srrz smz sn^z 



It may be mentioned that the derivative of p(z) is also an elliptic function 
with same periods and satisfy the relation 



z)=A\{[p{z)-e l }. (A.14) 



It is now straightforward to obtain equation (3) from equation (1) under the 
translation x = y^z + iK' by using the relations (A. 6) and (A.13). We will 
now apply the transformation y = (ex — p(z))/e 2 on equation (5). Noting the 
following relations 

p"{z) = 6p\z) - f , = 12p(z), (A.15) 

2 p {z) 

it is not very difficult to obtain the following equation 
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2 n [poo 



i=i 



dy 3 



+2^ [3 - m(m + l)]p(z) - 



'(l+l)e 2 e 3 
p(z) - ei 



e 2 



d 2 ^ 
cfo/ 2 



+ £ (e 2 ) 



+ < m(m + 1) 



l(l + l)e 2 e 3 
bW-ei] 2 



e 2 )^ = 



(A.16) 



The equation (7) will then readily follow under the transformation $(y) = 
[p(z) — ei] e ty on above equation by using the relations (A. 15). 
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